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Electrical  communications  at  room  temperature  have  a noise  background 
which  is  associated  with  thermal  (black  body)  radiation,  other  electromagnetic 
signals  such  as  lightning,  switching  transients,  and  electronics  noise. 

Nearly  all  geophysical  measurements  employ  an  electromagnetic  transducer. 

The  noise  is  partly  associated  with  the  amplifier-transducer,  partly  associated 
with  the  thermal  fluctuations  of  the  measuring  device,  and  partly  associated 
with  the  environment  background.  For  example  a gravimeter  is  employed  to  measure 
earth  surface  "g".  Such  an  instrument  is  usually  an  accelerometer  with  appropri- 
ate sensor  and  recorder.  Some  noise  originates  in  the  sensor  and  electronics, 
some  originates  in  the  accelerometer.  Usually  most  of  this  noise  has  its  origins 
in  local  earth  disturbances^  related  to  seismic  activity,  weather,  the  atmosphere, 
and  ocean  effects.  Gravity  gradiometers  have  sensitivity  limited  by  various  kinds 
of  electrical  and  mechanical  noise,  in  the  sensor  and  electronics. 

Noise  Limits 


If  there  were  no  significant  background,  and  the  electronics  were  noise  free, 
the  thermal  fluctuations  would  determine  the  ultimate  background  of  existing 
instruments.  As  temperatures  are  lowered  towards  absolute  zero,  the  thermal 
fluctuations  tend  to  zero.  Under  these  conditions  other  fluctuation  phenomena 
become  apparent.  These  are  the  zero  point  fluctuations  in  matter  and  certain 
quantum  noise  limits  in  electronics. 

The  ultimate  electronics  noise  arises  in  the  following  way.  Linear  elec- 
tronic amplifiers  employ  either  electrons  or  ions  in  excited  states.  These 
amplify  through  the  process  of  stimulated  emission.  In  the  absence  of  a signal, 

spontaneous  emission  may  occur.  This  is  a purely  random  process  and  contributes 
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noise.  The  temperature  corresponding  to  this  process  is  , 
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In  (1)  -tlw  Is  Che  energy  of  one  quantum,  and  k Is  Boltzmann's  constant. 
Intrinsic  Noise  of  An  Elastic  Solid 

If  a solid  is  in  thermal  equilibrium,  the  equlpartltlon  theorem  states 
that  each  normal  mode  of  oscillation  will  have  energy. 
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and  if  dLj  <<  1,  we  have 
kT  . 

d(w 


In  practice  ve  are  observing  strains  or  displacements  caused  by  some  effect 
which  we  are  attempting  to  measure,  and  the  strain  or  displacement  fluctuates 
with  an  amplitude  corresponding  to  the  normal  mode  energy.  The  extent  to  which 
these  thermal  fluctuations  interfere  with  measurements  depends  on  the  quality 
factor  Q,  in  the  following  way. 

For  angular  frequency  u>  there  is  a decay  time  t for  a transient  disturbance 
given  by 


r - (3) 

u 

The  energy  fluctuations  consist  of  noise  wavetrains  of  duration  roughly 
t seconds  before  the  amplitude  or  phase  changes  by  a large  amount,  on  average. 
If  we  make  some  observation  of  duration  At  the  thermal  energy  fluctuation  is 
expected  to  be^ 
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With  electronic  instrumentation  there  will  be  additional  noise  contributed 

by  the  amplifier  with  some  noise  temperature  I with  TM  > T„n  defined  by  (1)  . 

N N Ny 

The  minimum  "detectable  energy  change"  at  the  resonance  frequency  of  the  solid 
is  therefore 


| 

j 


kT  At 


kT, 


N 
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One  can  show  that  this  equation  holds  within  a factor  of  2 when  the  motion 

is  observed  below  its  resonance  frequency  with  At  2.  n/u.  The  quantities  Ac, 

» 

t,  and  T^  are  not  independent.  For  a given  strength  of  coupling  the  electronics 
to  the  system  being  observed  there  will  be  an  optimal  At.  Increasing  the  coupling 
may  decrease  the  Q and  thereby  decrease  T. 

Nonetheless  (5)  is  a good  guide  for  Indicating  how  the  technology  to  measure 
accelerations  and  gravity  gradients  is  improved  by  going  to  lower  temperatures, 
higher  quality  factor  materials,  and  lower  noise  electronics. 

All  presently  known  systems  are  described  very  well  by  (5).  As  temperatures 
are  lowered  further  and  quality  factors  also  improve,  eventually  (5)  becomes 
less  than  one  phonon  or  photon,  and  a number  of  modifications  appear  to  be 
required. 


As  already  noted,  if 


kT 


is  not  small,  kT  in  (5)  is  replaced  by 


kT 

e -1 


(5)  assumes  that  a system  is  employed  which  is  kept  in  thermal  equilibrium 
with  a heat  bath  at  temperature  T.  Historically  this  has  been  the  most  convenient 
way  to  operate  a system.  Also  the  linear  electronic  amplifiers  invented  and 
developed  in  this  century  have  all  been  stimulated  emission  devices  with  the 
limiting  noise  temperature  given  by  (1). 
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The  measurement  problem  can  be  set  up  in  different  ways  not  described  by 


these  equations.  For  example , the  quantum  mechanics  tells  us  that  certain 
quantities  can  in  principle  be  measured  with  arbitrarily  great  precision. 

Such  measurements  then  limit  the  precision  obtainable  with  measurement  of  other 
quantities  at  the  same  or  at  certain  later  times.  If  a system  is  in  thermal 
equilibrium  with  a heat  bath  only  its  temperature  is  known.  This  is  less  than 
the  information  required  to  define  a definite  quantum  state. 

Therefore  another  approach  to  the  measurement  problem  is  to  prepare  the 
system  so  that  it  is  in  a well  defined  quantum  state  rather  than  at  a definite 
temperature.  For  example,  one  might  drive  a system  with  a stable  oscillator 
aqd  detect  changes  by  observation  of  amplitude  or  phase  modulation  of  the  output. 
Quantum  mechanics  does  not  give  a clear  prescription  for  the  optimal  way  of 
preparing  a system  so  that  greatest  sensitivity  is  achieved,  and  discovery  of 
this  optimal  way  is  an  important  research  goal. 

Another  issue  is  the  electronics.  At  frequencies  from  d.c.  to  the  infra  red, 

stimulated  emission  devices  are  the  only  ones  available.  At  higher  energies 

particle  counters  may  be  employed.  A counter  may  be  in  its  ground  state,  emitting 

no  noise  whatsoever,  until  one  phonon  or  photon  appears  as  a result  of  some  change 
2 

in  a system.  No  successful  counter  of  this  type  has  ever  been  developed  for  the 
radio  region. 

We  may  summarize  these  comments  by  stating  that  systems  designed  in  the 
traditional  way  can  be  improved  by  going  to  lower  temperatures  and  employing  more 
nearly  perfect  crystals.  Eventually  the  energy  sensitivity  approaches  that  of  a 
single  photon  or  phonon.  In  principle  even  this  kind  of  performance  can  be 
exceeded  by  employment  of  systems  not  in  thermal  equilibrium  together  with  elec- 
tronics not  based  on  stimulated  emission. 
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Recent  Research 


By  1974  we  had  developed,  at  Maryland,  technology  involving  large  elastic 
solids  at  room  temperature.  T was  300  K,  At  was  1 second,  t was  40  seconds 
and  T^,  was  approximately  5 Kelvin.  We  could  therefore  observe  changes  in 
elastic  solid  strain  of  about  10  1 corresponding  to  energy  changes  of  about 
12  degrees. 

With  AFOSR  support  we  were  able  to  develope  a parametric^  capacitor  accel- 
erometer sensor  and  extend  our  technology  to  the  liquid  helium  region.  T was 
4.2  Kelvin,  the  fully  loaded  Q was  nearly  1,000,000,  and  we  could  observe 
changes  in  elastic  solid  energy  of  about  0.5  degrees.  Thus  we  improved  on  our 

i 

1974  technology  by  approximately  a factor  25. 

The  instrumentation  depends  in  an  important  way  on  the  quality  factor  Q 
of  the  elastic  solid.  During  the  past  year  we  studied  four  materials  — 

large  single  crystals  of  silicon,  sapphire,  lithium  niobate  and  ruby.  These 
results  to  be  described  later  illustrate  that  detailed  investigations  are  required 

to  realize  the  advantages  of  going  to  lower  temperatures,  for  practical  instrument 
designs . 

Cooling  a solid  removes  enormous  amounts  of  energy.  Sometimes  metastable 
states  are  "frozen  in"  and  a long  period  is  required  for  the  state  of  minimum 
free  energy  to  be  reached.  Excess  noise  and  poor  sensitivity  might  result  from 
the  release  of  stored  energy  at  a random  rate.  We  observed  such  effects  at 
liquid  nitrogen  temperatures,  but  not  at  liquid  helium  temperatures.  No  noise 
measurements  have  ever  been  made  on  the  very  high  Q single  crystals  at  liquid 
helium  temperatures  - by  very  high  Q we  mean  exceeding  Q values  of  ten  million. 
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An  attempt  was  made  Co  measure  the  noise  of  Che  large  silicon  crystal  at 
4 Kelvin  in  the  following  way.  The  crystal  is  relatively  soft  and  fragile. 

It  is  "dislocation  free".  For  these  reasons  we  did  not  wish  to  risk  permanent 
damage  by  use  of  epoxy  to  cement  a sensor  to  the  crystal  surface.  Therefore, 
a piezoelectric  crystal  was  bonded  to  an  aluminum  structure  which  served  as 
the  mount  for  this  crystal.  The  coupling  was  adequate  to  readily  excite  the 
silicon  crystal  by  driving  the  small  piezoelectric  crystal  with  a signal 
generator. 

A series  of  measurements  was  carried  out  during  the  past  few  months.  The 
quantity  which  determines  the  ratio  of  noise  from  the  elastic  solid  to  elec- 
tronics wideband  noise  is  8Q.  8 is  the  ratio  of  the  elastic  solid  energy 

appearing  at  the  output  terminals  to  the  internal  energy,  and  Q is  the  quality 
factor.  For  the  arrangement  of  Figure  2,  a 8Q  less  than  .3  was  observed,  with 
all  observed  noise  arising  entirely  in  the  electronics. 

We  have  a large  sapphire  crystal  which  is  extremely  rugged,  and  hard. 

Epoxy  and  other  bonding  materials  are  readily  removed  from  its  surface.  There- 
fore we  decided  to  obtain  a much  larger  8Q  by  bonding  a small  lead  zirconate 
titanate  crystal  to  the  large  sapphire  at  its  center.  This  was  successful  in 
raising  the  8Q  by  two  orders,  to  a value  5,  with  Q = 107.  This  value  of  8Q 
is  quite  adequate  to  carry  out  an  investigation  of  the  noise  originating  in  the 
sapphire.  Continuing  problems  with  our  liquifier  did  not  give  enough  stable 
operation  time  to  complete  these  measurements  before  March  1,  1979. 
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THE  MEASUREMENT  OF  NOISE  IN  MONOCRYSTALS  AT  CRYOGENIC  TEMPERATURES: 
EXPECTED  PROBLEMS  AND  PROPOSED  SOLUTIONS 
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April  19,  1979 
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Introduction 


It  is  planned  to  make  noise  measurements  of  monocrystals  of  sapphire 
and/or  silicon  near  the  thermodynamic  temperature  of  5 Kelvin.  This  report 
describes  some  of  the  expected  problems  of  making  such  noise  measurements 
with  high  elastic  Q materials  and  low  coupling.  Solutions  to  the  problems 
are  also  discussed. 

Table  I summarizes  the  properties  of  the  silicon  and  sapphire  crystals. 
Each  crystal  is  supported  by  a fourpoint  rigid  suspension  built  by 
Jean-Paul  Richard.  These  are  shown  in  figures* 1 and  2.  Each  is  coupled  to 
a small  piezoelectric  ceramic  transducer  (pzt-4) . For  the  silicon  the  trans- 
ducer is  expoxied  to  the  mount  as  shown  in  figure  1.  The  transducer  for  the 
sapphire  is  epoxied  directly  onto  the  crystal  itself. 

Figure  3 shows  the  well  known  equivalent  circuit  for  one  mode  of  an 
elastic  solid  coupled  to  a transducer  and  amplifier.  Resistors  R^  and  R£ 
couple  Brownian  motion  energy  into  the  crystal  mode.  This  energy  is  propor- 
tional to  the  thermal  temperature  T.  Other  processes  such  as  non-equilibrium 
effects  and  the  currert  noise  of  the  amplifier  also  put  energy  into  the 
crystal.  The  effective  noise  temperature  is  defined  for  the  purposes  of  this 
report  to  be  the  thermal  temperature  which  R^  and  R£  would  have  to  have  to 
produce  the  same  noise  energy  in  the  crystal  that  all  the  sources  combined 
produce.  It  is  the  goal  of  the  experiment  to  measure  the  effective  noise 
temperature  for  the  fundamental  mode  of  the  silicon  and/or  sapphire  crystals 
near  5 Kelvin. 


The  problems  in  making  the  noise  measurement  include  the  following.  High 
Q crystals  "ring"  for  a long  time.  The  noise  measurement  can  be  disrupted  if 

• : 

the  amplitude  of  the  ringing  greatly  exceeds  the  noise  level.  The  present 


crystal,  transducer,  and  amplifier  combination  requires  a very  small  bandwidth 
to  see  the  energy  noise  of  the  crystal  over  the  wideband  noise  of  the  ampli- 
fier. This  in  turn  requires  a very  stable  reference  frequency  tuned  accurately 
to  the  crystal  resonance.  A relatively  small  bandwidth  is  also  required  in 
one  of  the  intermediate  stages  of  amplification.  Before  making  the  noise 
measurements  electronics  must  be  carefully  designed,  built,  and  tuned  to  solve 
the  above  foreseeable  problems.  This  is  discussed  in  the  following  sections. 
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Summary  of  Noise  Analysis 


Consider  the  electronics  of  figure  4.  x(t)  and  y(t)  are  output  voltages 
which  are  slowly  varying  functions  of  time.  Let  G be  the  overall  gain  of  the 
electronics,  and  let  B(u>)  be  the  complex  transfer  function  of  the  Butterworth 
low  pass  filter  after  its  contribution  to  the  overall  gain  has  been  factored 


S(^)  = 


- 


B(o)  - 


/ + If  - CJ2^3 


4-  oo+'T* 


If  the  input  voltage  to  the  preamp  is  noise  only,  described  by  spectral  density 
ein(w).  then 
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where  wr  is  the  reference  frequency  of  the  lock-in  amplifier.  Further  linear 

filtering  of  x and  y,  either  analog  or  digital,  will  modify  the  above  noise 

equation.  For  filtering  relevant  to  this  report, 

r+ 


<(x+xr0)a+  = 


(^co  + U)r)  Jloo 


(3) 


It  can  be  shown  that  in  the  vicinity  of  the  resonance  of  the  crystal 
(figure  3)  the  voltage  spectral  density  at  the  preamp  input  is  given  by. 
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The  above  assumes  the  following: 


I co  I « CJ. 


(SO. 


COc  3 ^ a 


>>  / 


O'3/. 3 >>  / 


and  that  the  spectral  densities  of  the  noise  sources  are  slowly  varying  with 
frequency.  Notice  that  the  input  spectral  density  divides  naturally  into  two 
parts,  a wideband  part  (the  first  three  terms)  and  a narrowband  part  (term 
with  1 + u>2t02  in  denominator).  The  narrowband  term  is  proportional  to  the 
spectral  density  of  the  energy  of  the  crystal  for  the  mode  being  observed. 

The  noise  energy  of  the  crystal  comes  from  ej(ui),  e2(u0,  Jf(u>),  and  possibly 
other  sources  which  are  not  included  in  the  analysis.  If  thermal  Brownian 
motion  is  the  only  source  of  noise  then 


e3(^c)  + = JL  £t  (R,  +R*) 


The  effective  noise  temperature  (Tn)  of  the  Rj  + R2  combination  Is  defined  for 
this  report  to  be  the  thermodynamic  temperature  they  would  have  to  have  to 
produce  the  same  narrowband  noise  as  all  the  noise  sources  combined.  Thus, 


e3(oj „)  + e/(co.)  + 


CO?  cj 


(R'  + M 


0 + 

Spurt  es 


It  should  be  emphasized  that  this  noise  temperature  describes  only  the  narrow- 
band  part  of  the  noise,  those  processes  contributing  to  noise  energy  in  the 
crystal  mode  itself. 
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Let  e 2 (a»0)  be  the  sum  of  the  three  wideband  terms. 


C^o)  — ~h  C^o)  -h 


• a 


J f C^c) 


U>?  cf 


Then  equation  A for  the  input  spectral  density  can  be  rewritten, 

/ 


eZ,(UJ+OJ‘)  = ew(“M  + 


_ sjsaj,  7; 
/ ^ coar0a  ir  ooD  c a 


The  results  for  calculations  of  the  linear  filtering  given  by  equations 
1,  2,  and  3 are  shown  next.  It  is  assumed  that  the  reference  frequency  is 


equal  to  the  crystal  frequency,  u>r  » u>0. 
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No  approximations  were  made  in  evaluating  the  integrals  for  the  above. 


Note  the  following  relationship, 

<(x+  xr0)a-H  (y-  +■  y.7'.)*  > - 

<xa*ya>  + <xtpr.1 


Table  I gives  the  expected  values  for  the  elements  of  the 
equivalent  circuits  of  the  silicon  and  sapphire.  There  is  an  expected  error 
of  £ 20Z  for  the  sapphire  parameters.  Due  to  the  uncertainity  of  (5  for 
silicon  only  limits  exist  for  several  of  its  parameters.  The  voltage 
and  current  noise  of  the  FET  amplifier  to  be  used  has  not  yet  been  measured 
at  3468  and  18642  hertz.  At  1755  hz  the  noise  was  measured  as. 


efM  = 


_ 9 yol  + s 
l . 8 X / 0 


\farr 


— .72X10 


y/ol’t’S 


Jr-a-d.  / s e c 


j 


-IS  _/c 

•+-QXIO  _ 1.6  x 10  volts  — 

y 3.7T  slrc-l/sec. 

for  BF817  FET 


For  the  parameters  in  the  table  the  effect  of  the  current  noise  can  be  neglected. 
The  expected  noise  temperature  for  both  crystals  is  the  thermal  temperature  of 
4 to  5 Kelvin.  R2  has  not  been  directly  measured  yet  for  either  crystal.  If 
the  electronic  Q of  the  piezoelectric  transducers  is  estimated  conservatively 
at  50  then  the  effect  of  R2  on  the  Q of  the  crystal,  the  narrowband  noise,  and 

i k 

the  wideband  noise  can  be  neglected.  Next  the  time  constant  of  the  filters 


♦Mechanical  losses  in  the  piezoelectric,  the  epoxy,  and  the  suspension  contribute 
to  Rj . 
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^ will  be  chosen  so  that  Che  expected  wideband  and  narrowband  outputs  are 

equal.  Table  II  lists  the  values  for  tf  and  the  expected  total  noise.  These 
numbers  are  to  be  used  at  present  for  design  purposes  only.  After  better 
measurements  of  the  noise  and  parameters  are  made  there  could  be  considerable 
differences. 

To  obtain  the  narrowband  part  of  the  noise,  first  the  total  noise 

(wideband  plus  narrowband)  will  be  measured.  The  wideband  part  alone  can 

4* 

be  determined  by  slightly  detuning  the  reference  oscillator  of  the  lock-in 
amplifier  so  that  the  narrowband  part  falls  outside  the  overall  bandwidth. 
Substractlng  then  yields  the  average  narrowband  output.  From  this  measure- 
ment and  equations  5,  6 , 0r  7 the  noise  temperature  can  then  be  computed. 
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Crystal  Transient  Reduction 


During  a given  liquid  helium  run  the  crystal  will  be  excited  many  times 
above  the  noise  level  to  obtain  a measurement  of  its  frequency  and  Q.  An 
accurate  frequency  and  Q measurement  is  needed  for  noise  measurements.  How- 
ever, this  transient  can  interfere  with  the  noise  measurements.  Waiting  for 
it  to  die  away  could  take  a long  time  if  the  Q is  high.  The  amplitude  decay 
time,  tq,  is  8 1/6  hours  for  the  silicon  described  in  Table  I. 

For  the  sapphire  tq  - 990  seconds.  Thus  to  decay  by  a factor  of  106  the 
silicon  requires  4.7  days,  while  the  sapphire  needs  only  3.8  hours.  Therefore, 
measurements  of  the  silicon  noise  will  be  delayed  for  a substantial  amount  of 
time.  With  the  present  Q noise  measurements  of  the  sapphire  will  not  be 
substantially  delayed  by  its  excitation  transient.  However,  future  transducers, 
capacitive  or  inductive  instead  of  ceramic  piezoelectric,  will  by  design  produce 
much  higher  Q’s  for  all  crystals.  This  section  briefly  discusses  some  ways  to 
reduce  the  transient  itself  or  its  effects  on  the  noise  measurement. 

The  crystal  is  excited  by  applying  a sinusoid  voltage  at  the  crystal 
frequency  across  C2  (figure  3)  for  a short  interval  of  time.  The  transient 
can  then  be  reduced  by  applying  another  sinusoid  of  appropriate  phase  and 
duration.  This  technique  has  been  tried  successfully  with  the  silicon  bar  at 
liquid  helium  temperatures.  However,  it  does  not  enable  one  to  closely  approach 
the  noise  energy  of  the  crystal  due  to  the  relatively  large  amounts  of  wideband 
noise  compared  to  the  low  level  of  noise  energy  expected  in  the  crystal. 

Reducing  the  bandwidth  does  decrease  the  wideband  noise,  but  also  increases 
the  response  time  of  the  filters.  This  means  waiting  a long  time  to  get  infor- 
mation about  the  phase  and  amplitude  of  the  transient.  For  the  silicon 
described  by  Table  I as  the  transient  approaches  the  level  of  the  noise 
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energy  in  the  crystal,  the  time  constants  of  the  filters  needed  to  resolve  it 
exceed  the  damping  time  of  the  antenna.  Thus  this  method  will  not  eliminate 
the  lowest  levels  of  the  transient. 


A second  method  for  reducing  the  transient  is  to  increase  the  damping  by 
adding  an  external  resistor  across  C2.  However,  unless  C2  can  be  made  effec- 
tively smaller  (increased  coupling),  significant  energy  from  the  crystal  cannot 
be  dissipated  by  the  resistor.  This  can  be  done  by  adding  an  appropriate 
inductor  across  also.  In  order  to  reduce  the  transient  to  the  level  of  the 
energy  noise,  the  resistor  and  probably  also  the  inductor  will  have  to  be 
cooled  to  the  same  temperature  as  the  crystal.  Katsunobu  Oide  and  others  at 

the  University  of  Tokyo  have  used  feedback  to  reduce  the  noise  temperature  of 
1 

a resistor.  This  technique  might  be  applicable  to  removing  the  transient 
from  high  Q.low  coupling  crystals,  but  the  detailed  analysis  has  yet  to  be 
done  by  this  author.  After  the  transient  is  gone  the  external  resistor  and 
inductor  must  be  removed  to  make  the  noise  measurement. 

Electronic  filtering,  either  digital  or  analog  can  reduce  the  interference 
caused  by  the  transient  energy  in  the  crystal  when  noise  measurements  are  made. 
The  x and  y outputs  of  the  excited  crystal  have  the  following  form: 

- ±/T'o 

- -t/'T'o 


= 

(*) 

4- 

Kt  e 

11 

+ 

"n"  refers 

to  the  noise 

discussed 

the  subscript  "e"  refers  to  the  excitation  of  the  crystal  above  equilibrium 
at  some  initial  time.  The  ensemble  average  outputs  are  then  given  by. 


-.2  * /V 
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<■*'«)  + y.3U) > r/  = <X hV  £ > + (*/ i- ?: )5i 


3 a...  ‘M'C*  -a*/^ 

G 


<(xw+iwi'.)+  (?(*)+  ywn)5  > = 
<(x„  + x„r.)1+  (^»  + ^r.)a> 


Note  that  the  third  linear  filter  totally  removes  the  excitation  transient  of 
the  crystal  before  squaring.  In  practice  the  crystal  decay  time,  t0,  will  not 
be  known  well  enough  to  remove  all  the  transient.  Suppose  r ' is  the  value 
which  is  used  for  tq.  Then  the  following  term  must  be  added  to  the  latter 
equation  above,  „ . 0 

(*•*<-  ?.■)(,  -2Lf 

For  the  present  silicon  and  transducer  one  or  more  of  the  above  methods 
can  be  used  to  remove  the  transient  if  an  attempt  at  measuring  the  energy 
noise  is  made.  Increased  coupling  will  improve  the  effectiveness  of  these 
techniques. 
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Reference  Frequency 


The  lock  in  amplifier  requires  a reference  oscillator.  The  results  of 
the  noise  calculations  — equations  5,6,  and  7 — are  valid  only  for  the 
reference  frequency  sufficiently  close  to  the  crystal  resonance.  When  the 
filter  time  constant  is  roughly  equal  to  the  decay  time  of  the  crystal,  it 


turns  out  that  for  all  three  of  the  linear  filters  considered  the  following 


condition  must  be  maintained. 


The  crystal  frequency  must  be  measured  to  better  than  this  accuracy.  The 
reference  must  be  tuned  to  and  maintained  at  this  frequency  to  better  than 
this  accuracy  throughout  the  duration  of  the  measurement.  Data  must  be 
collected  for  a time  greater  than  the  decay  time  of  the  crystal.  Thus  the 
stability  and  accuracy  requirements  for  each  of  the  crystals  is. 


silicon,  accuracy  and  stability  of  better  than  one  part 
in  3 x 10®  over  a time  greater  than  8 hours. 

sapphire,  accuracy  and  stability  of  better  than  one  part 

•j 

in  6 x 10  over  a time  greater  than  hour. 

Devices  are  commercially  available  which  have  the  needed  stability  and 
accuracy.  They  cost  several  thousands  of  dollars  to  buy  and  several  hundred 
of  dollars  per  month  to  i nt.  These  alternatives  might  be  considered. 

Frequency  measuring  and  generating  devices  of  sufficient  accuracy  and 
stability  have  been  constructed  from  equipment  which  we  already  have  in  our 
possession.  The  Hewlett  Packard  10544A  lOMhz  quartz  crystal  oscillator  has 
a rated  stability  of  better  than  one  part  in  2 x 109  per  day.  The  oscillator 
and  the  HP  5216A  electronic  counter  have  the  needed  precision  to  measure  the 
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crystal  frequencies.  The  lOMhz  oscillator  is  also  used  to  generate  the 
reference  frequency.  This  eliminates  the  need  for  calibration.  Figure  5 
shows  the  phase  lock  loop  which  has  been  used.  Note  that  the  lOMhz  oscillator 
is  not  in  the  loop.  A second  crystal  oscillator  is  servoed  by  a voltage  which 
is  proportional  to  the  difference  between  the  desired  and  the  generated  periods. 
Prototypes  of  this  frequency  synthesizer  have  been  successfully  tested  for  both 
silicon  and  sapphire  frequencies. 

The  frequency  of  the  crystal  is  a function  of  temperature.  McGuigan  and 

others  have  determined  that  at  liquid  helium  temperatures  the  temperature 

2 

dependence  of  the  resonance  frequency  of  silicon  has  the  following  form, 

~ ^o)  — i.  8 a x icT'°t3 
■fCo) 

During  a typical  run  with  the  cryostat  to  be  used  to  make  our  noise  measurements 
the  temperature  of  the  experiment  changed  from  4.4  Kelvin  to  5.7  Kelvin  over 
the  six  day  holding  period  of  the  dewar.  The  best  temperature  stability  was 
about  .1  Kelvin  per  day  for  the  first  four  days.  An  adequate  noise  measurement 
requires  collecting  data  for  a time  equal  to  many  decay  times.  Ten  decay  times 
for  the  silicon  is  3.4  days.  During  this  time  the  frequency  of  the  silicon 
will  drift  by  about  one  part  in  3 x 108.  Thus  the  best  temperature  stability 
is  Just  at  the  edge  of  what  is  needed.  Improvements  have  recently  been  made  in 
the  cryostat  which  are  expected  to  increase  the  thermal  stability,  ie.  better 
heat  sinking  of  the  wires  going  to  the  experiment.  However,  this  has  not  yet 
been  tested.  Active  temperature  control  can  be  provided  if  necessary.  Alter- 
natively if  the  temperature  dependence  of  the  frequency  is  known  accurately 
enough,  the  reference  frequency  can  be  compensated  to  follow  the  crystal 
frequency. 
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The  temperature  dependence  of  the  resonance  frequency  of  sapphire  has 
not  yet  been  obtained  by  this  author.  Such  data  has  probably  been  obtained 
by  V.  Braginsky.  We  hope  to  obtain  it  in  the  near  future  from  him  or  our 
own  measurements.  However,  the  frequency  stability  requirements  for  the 
sapphire  presently  being  tested  are  less  stringent  than  for  silicon  by  a 
factor  of  over  100.  Thus  noise  measurements  for  the  sapphire  probably  won't 
need  improved  temperature  stability  or  compensation. 
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Cain  and  Bandwidth 


Careful  choice  of  gain  and  bandwidth  of  the  stages  of  the  electronics 
of  figure  4 is  needed  to  obtain  an  adequate  output  without  saturation. 
Consider  first  a specific  example,  x2  + y2  for  the  sapphire.  For  a .7  volt 
rms  output  for  each  of  x and  y a total  gain  of  1010  is  needed.  The  gain 
in  the  low  pass  filters  should  be  no  more  than  100  due  to  dc  drifts.  Thus 
the  gain  of  the  preceeding  stages  should  be  greater  than  108.  The  detectors 
of  the  lock  in  amplifier  (a  PAR  129)  begin  to  saturate  at  about  one  volt  rms 
of  noise.  The  narrowband  signal  will  not  cause  these  detectors  to  saturate 
but  the  wideband  noise  might.  The  effective  noise  bandwidth  of  the  preamp 
and  postamp  which  will  not  produce  this  saturation  is  given  by 

io8*  s x < 1 i/olf 

assuming  2 nv//  hz  wideband  noise  at  the  preamp  input.  Thus  we  need 
Af  < 25  hertz.  This  narrow  bandwidth  at  18642  hertz  might  just  barely  be 
attainable  by  the  four  cascaded  LC  tuned  amplifiers  of  the  postamp.  This 
device  presently  has  a bandwidth  of  about  78  hertz.  More  careful  tuning  and 
slight  amounts  of  positive  feedback  might  bring  the  bandwidth  down  some  more. 

Table  III  gives  the  expected  upper  limit  for  the  postamp  bandwidth  for 
silicon  and  sapphire  and  each  of  the  three  linear  filters.  With  the  possible 
exception  of  sapphire  at  25  hz  = Af,  amplifiers  with  quartz  crystal  filters 
or  a combination  of  down  conversion  and  crystal  filters  will  probably  be 
needed  to  achieve  these  narrow  bandwidths. 
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Conclusion 


The  noise  measurement  for  the  sapphire  appears  feasible  with  a moderate 
effort  put  into  further  electronic  development.  The  major  work  which  remains 
to  be  dor.e  is  the  development  of  a postamp  with  narrow  enough  bandwidth.  It 
is  well  within  technical  feasibility. 

The  noise  measurement  for  the  silicon  is  more  difficult,  due  mainly  to 
the  low  coupling.  Much  more  development  is  needed  than  for  the  sapphire. 

It  might  be  more  feasible  to  improve  the  coupling  instead.  This  would  relax 
the  present  requirements  for  transient  removal,  frequency  and  temperature 
stability,  and  bandwidth  by  giving  a larger  energy  noise  signal  and  permitting 
a smaller  filter  time  constant  to  be  used. 

Larger  couplings  can  be  achieved  with  piezoelectric  transducers  and 
epoxy,  but  this  reduces  the  Q.  As  a consequence  the  Brownian  motion  noise 
of  the  crystal  comes  more  from  loss  processes  in  the  transducer  than  from  the 
crystal  (and  its  suspension).  The  sapphire  presently  being  used  had  a Q of 
A x 10B  at  4 Kelvin  before  the  piezoelectric  was  attached.  The  present  Q 
of  58  x 106  means  that  at  equilibrium  roughly  15X  of  the  energy  noise  to  be 
observed  comes  from  the  crystal  (and  its  suspension).  Possible  future 
experiments  with  these  crystals  could  have  8's  as  large  as  10-7  with  inductive 
and  capacitive  transducers  while  maintaining  high  Q.  These  would,  however, 
require  separations  of  about  .1  mm  between  the  transducer  and  the  end  of  the 
crystal.  The  use  of  a properly  matched  SQUID  amplifier  instead  of  an  FET 
would  also  be  a significant  Improvement. 
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PROPERTIES  OF  MONOCRYSTALS  COUPLED  TO  TRANSDUCERS 
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Figure  3,  Equivalent  Circuit  and  Definitions 
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Figure  4,  Electronics 

x and  y may  be  recorded  directly  on  magnetic  tape  or  undergo 
further  analog  processing  before  being  recorded. 
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Q MEASUREMENTS  AT  CRYOGENIC  TEMPERATURES 


As  noted  in  equation  (5),  the  ultimate  sensitivity  which  can  be  attained 
in  measurements  of  changes  in  energy  of  any  elastic  solid  detection  depends 
on  the  free  oscillation  Q --decay  time  t.  A program  to  study  the  low  tempera- 
ture noise  must  begin  with  Q measurements.  The  details  of  the  mount,  and 
instrumentation  will  depend  critically  on  the  Q values. 

During  the  past  year  detailed  measurements  of  the  Q of  the  fundamental 
longitudinal  mode  of  ruby,  sapphire,  lithium  niobate  and  silicon  were  carried 
out.  All  crystals  were  cylindrical  rods.  Most  of  the  measurements  were  made 
supporting  the  crystals  by  aluminum  four  point  suspensions.  One  measurement 
(ruby)  was  made  with  a tungsten  wire  suspension  as  shown  in  Figure  1.  Either 
end  capacitor  plate  or  piezoelectric  ceramic  crystals  were  employed. 

Ruby  is  sapphire  doped  with  (paramagnetic)  chromium  ions.  It  is  being 
studied  in  order  to  assess  the  importance  of  paramagnetic  ions  in  affecting 
the  Q.  Most  materials  contain  a certain  amount  of  paramagnetic  ions  as 
impurities.  At  sufficiently  low  temperatures  the  specific  heat  of  the  lattice 
vibrations  of  a crystal  will  be  extremely  small  — so  small  that  measurements 
would  be  affected  by  the  heating  effect  of  the  high  energy  cosmic  rays. 

Control  of  the  specific  heat  by  doping  with  paramagnetic  impurities  may  be 
a very  useful  technique. 

The  instrumentation  of  a nearly  perfect  crystal  to  observe  its  low  temp- 
erature properties  is  enormously  simplified  if  the  crystal  is  piezoelectric. 

It  is  therefore  of  interest  to  explore  single  crystal  materials  such  as 
lithium  niobate. 
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RUBY 


Two  differenc  ruby  crystals  have  been  tested.  Each  is  a single  crystal 
of  pink  sapphire,  sapphire  doped  vith  . 05Z  chromium.  They  both  have  a rough 
polish  on  the  sides  and  ends.  Their  mass  and  dimensions  are; 

.96  kg  mass,  3.0  cm  diameter,  33.1  cm  length 
1.44  kg  mass,  5.0  cm  diameter,  18.5  cm  length 
The  .96  kg  ruby  was  suspended  with  a tungsten  wire  as  shown  in  figure  1. 
A pzt-4  piezoelectric  ceramic  was  bonded  to  the  surface  of  the  crystal  with 
Stycest  2850  FT.  The  Q of  the  fundamental  longitudinal  mode  at  4.2  Kelvin 
was  20  x 106.  Figure  5 shows  Q measurements  between  30  and  300  Kelvin. 

• The  1.44  kg  ruby  was  supported  vith  the  four  point  aluminum  suspension 
shown  in  fij-ure  2.  The  end  surfaces  were  coated  with  aluminum,  so  that  a 
split  capacitor  plate  could  be  used  to  measure  the  motion.  This  is  shown  * 
in  figure  3.  The  Q of  the  fundamental  longitudinal  mode  at  4.3  Kelvin  was 
47  x 106. 

Results  reported  in  the  next  sections  suggest  that  the  present  aluminum 
four  point  suspension  reduces  the  Q more  than  a wire  suspension.  If  this  Is 
true  then  the  lover  Q of  the  wire  supported  ruby  would  then  be  due  to  the 
piezoelectric  transducer. 
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SAPPHIRE 


The  sapphire  crystal  vas  obtained  from  Union  Carbide.  It  is  25  cm 

long  and  8. 25. cm  in  diameter  and  has  a mass  of  5.2  kg.  The  side  surface 

vas  unpolished  and  untreated  in  any  way.  The  end  faces  had  a rough  polish 

and  were  coated  with  aluminum.  The  four  point  suspension  and  capacitor 

transducer  vere  the  same  as  that  used  for  the  1.44  kg  ruby,  figures  2 and 

3.  The  Q cf  the  fundamental  longitudinal  mode  at  4.3  Kelvin  vas  about 

400  x 10*.  Douglass  has  reported  measurements  of  the  Q of  the  quadrapole 

mode  of  a 5.6  kg  sapphire  disk  at  4.3  K.^  The  crystal  vas  supported  by 

two  ball  bearings,  centered  top  and  bottom.  The  result  vas  also 

Q - 400  x 10*.  Braginsky  has  also  published  measurements  of  the  Q of  a 
2 

sapphire  crystal.  This  crystal  vas  a cylindrical  rod  supported  at  its 
center  vith  the  wire  suspension  of  figure  1.  The  result  for  the  longi- 
tudinal mode  at  4.3  Kelvin  vas  5 x 109.  The  vire  suspension  produced  a 
Q vhich  vas  an  order  of  magnitude  higher  than  the  tvo  rigid  suspensions. 
Factors  such  as  crystal  purity  and  the  transducer  as  veil  as  suspension 
could  have  lowered  the  Q.  If,  however,  mechanical  losses  were  limited  by 
the  rigid  suspension,  then  the  use  of  higher  Q materials,  such  as  sapphire 
itself,  might  reduce  these  losses.  The  aluminum  alloy  which  has  been  used 
at  Maryland  for  all  rigid  suspensions  is  alloy  6061.  This  alloy  has  a 
mechanical  Q of  about  106  at  liquid  helium  temperatures.  Hirakawa  has 
discovered  that  aluminum  alloy  5056  has  a Q of  70  x 106  belov  20  Kelvin. ^ 
Our  future  research  into  suspensions  for  sapphire  crystals  will  Include 
four  point  suspensions  made  of  high  Q materials  such  as  aluminum  5056  and 
sapphire. 
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Tn*  Q's  of  the  5.2  kg  sapphire  and  the  1.44  kg  ruby  were  measured 
vi:h  the  *it«  transducer  and  suspension.  The  result  for  ruby  was  8>i  times 
lever.  This  suggests  significant  attenuation  due  the  chromium  ions  in 
the  ruby  material. 

I. 


SILICON 


Tvc  different  silicon  crystals  were  tested.  The  smaller  of  the  two 
was  furnished  by  Texas  Instruments.  Its  mass  is  .78  kg,  diameter  5 cm, 
and  length  17  cm.  The  sides  were  etched  and  the  ends  unpolished.  The  Q 
of  the  bending  code  at  11337  hertz  vas  observed  by  mounting  a pzt-4  piezo- 
electric transducer  under  one  of  the  four  aluminum  suspension  points  shown 
in  figure  2.  The  Q as  a function  of  temperature  is  shown  In  flugre  6.  At 
4.3  Kelvin  the  Q was  76  x 106. 

Data  fer  the  15  kg  silicon  crystal  furnished  by  Monsanto  are  given  on 

the  next  page.  The  surface  of  the  crystal  has  not  yet  been  improved.  The 
« 

ends  and  the  sides  have  visible  irregularities  resulting  from  the  crystal- 
lization process.  The  crystal  vas  supported  by  the  aluminum  four  point 
suspension  shown  in  figure  4.  A pzt-4  piezoelectric  ceramic  vas  epoxied 
onto  the  suspension  as  shown  in  this  figure.  Figure  7 shows  the  results 
of  Q measurements  as  a function  of  temperature^  At  5.8  Kelvin  the  Q was 
320  x 106. 

Both  the  above  silicon  crystals  gave  a significant  drop  in  Q in  the 
10  to  15  Kelvin  range.  The  15  kg  silicon  shows  a minimum  at  12  Kelvin. 

A minimum  ir.  the  Q of  silicon  (at  13  Kelvin)  has  also  been  observed  by 

4 

McGuiga  et  al  of  the  University  of  Rochester.  Their  4.9  kg  silicon  crystal 
was  polished  and  coated  with  aluminum  on  the  ends.  A split  capacitor  plate 
vAs  used  at  the  transducer.  The  sides  were  not  altered.  The  crystal  vas 
suspended  with  a tungsten  wire  and  Q measurements  were  made  from  1.3  to 
300  Kelvin.  The  minimum  in  the  Q curve  for  two  different  silicon  crystals, 
with  different  suspensions,  transducers,  and  end  face  conditions  strongly 
indicates  that  the  mechanism  causing  the  friction  comes  from  a bulk  property 
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SILICON  CRYSTAL 


15  kg  SILICON  CRYSTAL 


Crystal  Data  *99-0491 


Weight : 

Ci ameter : 

Material : 

Length : 

Crr.ductivity  Type 
F-es  istivitv : 

Major  Ir.puri.ties: 


Method  of  Growth: 


15,465  grams 

78.4  mm 

Single  Crystal  Silicon 

53.5  inch 
P-type 

20  ohm  cm  (seed  end)  - 14  ohm  cr.  (opposite 

Oxygen:  - 10le  atoms/cc 

Carbon:  - 1C17  atoms/cc 

Boron:  ~ S x 10  V*  atoms/cc 

Czochralski  (pulled  from  melt  contained  in 
quart:  crucible) 


Rate  of  Crystallization: 

100  mm/hr. 


This  very  large  crystal  is  described  as  "dislocation  free". 
Industry  employs  such  a description  for  crystals  having  less 
than  10  dislocations  per  square  centimeter. 


en 
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of  Che  silicon.  Observations  of  Che  silicon  bending  mode  data,  although 
they  don't  shew  a minimum  due  to  insufficient  data,  also  support  this 
conclusion.  Both  the  Rochester  4.9  kg  and  Maryland  IS  kg  crystals  were 
supplied  by  Monsanto.  Both  are  p-type  boron  doped  with  a similar  level 
of  other  impurities.  The  Rochester  group  suggests  that  it  is  the  boron 
impurity  which  produces  the  Q minimum.  It  would  be  most  interesting  to 
make  similar  measurements  for  a n-type  silicon  crystal. 

At  5.8  Kelvin  the  Maryland  15  kg  silicon  crystal  had  a Q of 
320  x 10s.  At  the  same  temperature  the  Rochester  group  reports  a Q 
over  three  times  higher.  This  difference  is  probably  due  to  losses 
from  the  suspension  and/or  piezoelectric  transducer.  Four  point  sus- 
pensions of  higher  Q will  be  tried  in  future  experiments.  Possible 
materials  include  aluminum  5056,  silicon,  and  quartz.  The  latter  two 
materials  offer  the  possibility  of  an  increased  number  of  chemical  bonds 
between  the  surfaces  of  the  crystal  and  the  suspension,  and  presumably 
lower  losses  at  these  points.  The  sides  and  ends  of  the  15  kg  crystal 
will  be  ground  and  polished,  and  a capacitor  transducer  installed. 

There  is  one  point  in  the  Q curve  for  the  15  kg  silicon  which  is 
significantly  lover  than  then  the  others  near  the  same  temperature 
(16  Kelvin  of  figure  7) . It  is  possible  that  this  Q was  measured  with 
an  amplitude  much  greater  than  the  other  measurements.  Alternatively 
it- may  have  a bad  measurement.  The  data  is  now  being  carefully  studied. 


CONCLUSION 


Many  Q measurements  of’ monocrystal*  of  ruby,  sapphire,  and  silicon 
at  cryogenic  temperatures  are  presented.  Comparison  of  these  results 
with  that  of  other  groups  suggests  that  higher  Q's  can  be  achieved  with 
rigid  suspensions  of  higher  Q material. 

The  minimum  in  the  Q curve  of  silicon  between  10  and  15  Kelvin  first 
observed  by  the  group  at  the  University  of  Rochester  has  been  duplicated 
in  a crystal  of  similar  composition  but  different  suspension,  transducer, 
and  end  surfaces.  This  strongly  supports  the  Rochester  group's  suggestion 
that  mechanical  losses  in  silicon  at  these  temperatures  are  due  to  internal 

t 

bulk  properties  of  the  material. 
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.78  kf»  silicon  crys 


LITHIUM  NIOBATE  AT  CRYOGENIC  TEMPERATURES 

•The  planar  piezoelectric  coupling  constant  of  a small  lithium  niobate 
disc  was  measured  at  300  K,  77  K,  and  4.2  R.  The  result  was  about  10"S  for 
all  three  temperatures.  This  agrees  with  a theoretical  calculation  using 
th*  piezoelectric  and  elastic  constants  of  lithium  niobate.  The  frequency 
v**  sbout  200  khz.  The  overall  Q vas  also  measured  at  4.2  K and  found  to 
be  23,000  to  30,000. 
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The  lithium  niobate  (LiNbO^)  crystal  was  obtained  from  Crystal  Technology 
Inc.  It  vas  a disc  with  diameter  .79  Inches  and  thickness  .05  Inches.  Silver 
vas  vacuum  deposited  on  each  side.  A small  area  of  the  silver,  about  two 
square  millimeters,,  vas  removed  from  the' center  of  each  side  of  the  disc.  Two 
AV3  26  tinned  copper  vires  were  attached  with  conducting  epoxy  TRA-CON  BB-2902. 
The  epoxy  made  a bond  with  the  crystal  and  with  the  silver.  Figure  1 shows 
hcv  the  crystal  vas  mounted  in  a small  vacuum  chamber.  The  chamber  vas 
^rmersec  in  liquid  nitrogen  or  liquid  helium  in  a small  glass  devar. 

Figure  2 shows  the  circuit  used  to  measure  the  electromechanical  coupling 
of  the  crystal.  The  capacitance  of  the  coaxial  line  from  A to  B vas  measured 
with  the  crystal  and  the  capacitance  box  disconnected. 

Figure  3 shows  the  equivalent  circuit.  Only  two  sources  of  loss  are 
shown  in  the  circuit.  There  are  certainly  others.  The  cables  going  from 
the  oscillator  to  signal  in  and  from  signal  out  to  the  oscilloscope  each  have 

« ft  N4 

a capacitance  of  about  100  pf.  Hence,  from  point  A the  signal  in  and  signal 
outdare  essentially  grounded.  Let  Z be  the  impedance  of  point  A to  ground. 

The  amplitude  of  the  signal  observed  on  the  oscilloscope  will  be  proportional 
to  jZ.J  • For  the  high  Q approximation,  the  signal  will  be  a minimum 
at  .frequency. 
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The  observed  signal  vill  be  a maximum  at  frequency. 


cOa  — 


CO 


m tf-jr 


L, 


Tne  electromechanical  coupling  is  defined  as 

C,  ' 

P ~ c.x 

Tne  inverse  of  this  is  sometimes  called  tbe  capacitance  ratio  or  the  piezo- 
electric ratio.*  can  be  calculated  from  the  above  equations  if 

^ j Cc  J Ca/nd  C3  are  toown-  The  result  is. 


P = 


ace  ■+-  c-»  +c 


CJ 


rHCLJC 


, ^ a 

>m'h 


— I 


£<5  • and  „ were  measured  with  the  frequency  counter  at  the  minimum 

and  maximum  signals.  was  measured  in  the  following  way.  When  the 

oscillator  was  set  at  a frequency  away  from  resonance,  the  impedance  7 
was  due  essentially  to  SCc  + Cj-h  C*  . The  signal  amplitude  on  the 
oscilloscope  was  recorded.  The  capacitance  box  was  removed  from  the  coaxial 
line  at  A,  figure  .2.  A capacitor  was  attached  to  the  box  in  place  of  the 
coax  and  found  to  give  the  same  signal  amplitude.  The  capacitance  of  the 
capacitor  and  its  connector  was  measured  with  a bridge.  It  was  80  pf. 
Tnerefore  Cs  -h.Cj  " 80pf,  Ca“  pf. 


* Cady,  Walter  G.,  Piezoelectricity,  p 353. 
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T-s  resonances  of  the  crystal  were  observed  from  100  khz  to  5 Khz. 
ite  res2r2r.ce  near  200  khz  had  changes  of  Impedance  much  larger  than  any 
tr.e  other,  resonances.  The  table  shovs  the  results  of  measurements  and 
calculations  for  this  resonance.  At  liquid  nitrogen  and  liquid  helium 
te=Pftr«--ures  y3  ■ 10  3.  Rt  and  /?a  can  be  estimated  from  this  data. 

At  CO  — ’ i*  m 


2L  ~ R-, 


70SI 


OJ  ■=  Ui 
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z: 


m + K 


-h 
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R?  + 


=> 


S'tfO.Sl 


The  Q of  this  mode  of  the  crystal  and  its  mount  is  then 


a = 


1 


coc, 


*+•  ^7* 


£ Sj  OOO 


The  Q cf  the  resonance  was  also  measured  at  liquid  helium  temperature 
by  sweeping  the  oscillator  quickly  through  the  resonance  and  recording  the 
transient  response  with  a memory  oscilloscope.  The  result  vas  about  30,000. 


i.ie  previous  experimental  results  will  now  be  compared  to  theory. 

__  * 

The  planar  coupling  constant  for  a piezoelectric  disc  is  given  by 

» » 

2 

2 

P ■ l:p  • (rationalized  MRS) 


p - k?“ 

m — 

cT 

C33 

to 

♦ « 

From  room  temperature 

data  on  ! 

LiKbOj 

Incorporated, 

d31 

.085 

X 

Hf11 

C/N 

T 

c-*3  " 

.254 

X 

IQ’9 

F/m 

r 

5.83 

X 

io-12 

m2/N 

S12  " - 

■1.15 

X 

io-12 

m2/N 

Thus,  kp2 

- 8 

m 

1.22 

x 10"' 

This  agrees  with  the  room  temperature  measurement  of  1.2  x 10  for  the 
coupling.  At  liquid  nitrogen  and  helium  temperatures  the  coupling  was 
1.0  x 10  3 . It  is  not  clear  if  this  difference  is  due  to  error  in  the 
technique  or  real  change  in  the  coupling.  However,  for  most  practical  purposes 
one  may  conclude  that  the  coupling  does  not  change  significantly  from  room 
temperature  down  to  cryogenic  temperature. 


IRE  Standards  on  Piezoelectric  Crystals:  Measurements  of  Piezoelectric 
Ceramics,  1961.  This  standard  is  reprinted  in  the  following  book. 
Piezoelectric  Ceramics,  Jaffa,  B.,  Cook,  V.  C. , and  Jaffe,  H.,  Academic 
Press,  New  York,  1971,  See  page  291  for  equation  for  kp2. 
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The  previous  planar  coupling  constant  is  for  the  radial  mode  of  the 

^-sc*  "Ts  ensure  that  this  mode  was  actually  observed,  the  expected  resonant 

• • 

r * 

frequency  vas  calculated  from  the  known  properties  of  the  lithium  niobate 
crystal.  The  result  is  208  khz  which  agrees  reasonably  with  the  observed 
frequency  of  196  khz  at  room  temperature.' 


Jaffe,  Cook,  and  Jaffe,  p. 


32 
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LARGE  20  nK  CRYOSTAT 


Figures  8,  9,  10,  and, 11  show  a cryostat  which  is  designed  to  cool  a 
15  kg  silicon  crystal  to  20  mK  and  maintain  it  there. 

The  cryostat  is  presently  under  construction  at  the  Gravity  Research 
Building  located  just  off  the  University  of  Maryland  campus  in  a relatively 
quiet  and  isolated  environment.  This  laboratory  already  contains  a working 
closed  cycle  liquid  helium  system,  consisting  of  a CTI  model  1400  liquifier, 
associated  compressors,  sod  a helium  storage  tank.  . The  system  contains  enough 
gas  to  fill  an  appropriately  sized  dewar  with  about  200  liters  of  liquid  helium. 

this  system  we  can  produce  a large  4.2  Kelvin  environment.  To  produce 
cillikelviz.  temperatures  we  are  presently  modifying  a large  vacuum  chamber  and 
dewar  to  accept  a SHE  DRI-423  helium  dilution  refrigerator  which  is  expected 
to  be  delivered  early  in  1979. 

A description  of  the  present  design  of  the  cryostat  will  now  be  given. 

The  cryostat  is  suspended  Inside  a large  vacuum  chamber,  10  feet  long  by  8 feet 
in  diameter.  The  copper  vapor  shield  hangs  from  the  inside  top  of  the  chamber 
through  acoustic  Isolation  filters.  Boil  off  liquid  helium  vapor  from  the  dewar 
circulates  through  copper  tubing  soldered  to  the  vapor  shield,  thus  cooling  it. 
This  shield  is  covered  with  several  layers  of  aluminum  coated  mylar  and  nylon 
net  (super insulation)  for  radiation  shielding.  The  liquid  helium  dewar  is 
enclosed  by  the  vapor  shield,  but  hangs  from  its  own  acoustic  isolation  filters 
attached  to  the  vacuum  chamber,  as  shown  inthe  figures.  The  stainless  steel 
dewar  is  connected  to  the  Isolation  via  a dacron  strap.  The  thermal  conduction 
down  this  strap  is  expected  to  be  the  main  heat  leak  into  the  dewar.1  If  the 
upper  end  of  the  strap  is  thermally  shorted  to  the  vapor  shield  then  the  expected 
boil  off  rate  is  12  liters  per  day  out  of  a total  of  200  liters  in  the  dewar. 

The  devar*  s outside  and  inside  surfaces  are  covered  with  1/16  inch  copper  shells 
to  provide  uniform  temperature.  Covering  the  outer  copper  shell  is  several 
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layers  of  superinsulation.  Inside  the  devar  is  a copper  shield  to  be  maintained 
at  1.3  Kelvin  by  the  cold  plate  of  the  dilution  refrigerator.  Rods  made  of 
nylon  or  vespel  SP-22  (Dupont  graphite  filled  polyimide  resin)  provide  the 
necessary  support  of  the  1.3  Kelvin  shield  and  thermal  isolation  from  the  4.2 

i % 

Kelvin  copper  shell  .of  the  devar.  Inside  the  1.3  Kelvin  shield  are  tracks 
upon  which  a cart  supporting  the  rest  of  the  experiment  can  be  rolled  in  and 
out.  As  shown  in  figures  9 and  10,  acoustic  isolation  filters  (alternating 
layers  of  felt  and  steel)  sit  on  the  cart.  These  filters  support  a platform 
which  is  also  maintained  at  1.3  Kelvin.  The  15  kg  silicon  crystal  is  coated 
with  niobium  and  supported  inside  a 20  mK  copper  shield  by  magnetic  levitation. 
The  latter  requires  superconducting  colls  with  persistent  currents  mounted 

i 

inside  the  20  mK  shield.  This  shield  rests  on  a ring  suspension  which  is 
supported  on  the  1.3  K platform  by  vespel  rods.  A strong  thermal  link  is 
made  from  the  ring  suspension  to  the  mixing  chamber  of  the  dilution  refriger- 
ator with  an  extremely  pure  copper  braid.  The  expected  heat  flow  from  the 
1.3  K platform  to  the  20  mK  ring  suspension  through  the  vespel  rods  is  about 
3 nicrovatts.  Other  expected  sources  of  heat  leak  are  less  than  microwatt. 
The  dilution  refrigerator  can  maintain  the  experiment  at  20  mK  with  a heat 
leak  of  3%  microwatts  with  a 2 mK  temperature  drop  along  the  copper  braid. 

This  assumes  that  the  heat  generated  by  the  persistent  current  in  the  coils, 
eg.,  from  flux  jumps,  is  smaller  than  the  heat  leak.  Further  design  and 
research  might  be  needed  to  ensure  this. 

There  are  two  separate  vacuums  in  the  cryostat.  The  space  between  the 
vacuum  chamber,  the  vapor  shield,  and  the  devar  is  kept  evacuated  with  a 
diffusion  pump,  this  space  contains  the  superinsulation  which  needs  a good 
vacuum  to  be  effective.  The  4.2  Kelvin  copper  shell  Inside  the  devar  is  the 

e« 

second  vacuum  space.  It  has  its  own  pumpout  line.  The  1.3  Kelvin  and  20 
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mil  ikelvin  shields  are  not  vacuum  ClghC.  During  cooldown  several  microns  of 
pressure  of  helium  will  be  ‘left  In  this  space  to  allow  cooling  of  everything 
by  conduction.  At  4.2  Kelvin  this  gas  is  pumped  out  and  further  cooling  Is 
accomplished  by  the  dilution  refrigerator.  The  levitation  colls  are  not 
turned  on  until  the  crystal  has  been  cooled  to  20  mK.  The  crystal  has  to  bs 
levitated  slowly  to  enstira  that  the  changing  magnetic  field  does  not  heat  It 
up. 

The  dilution  refrigerator  must  be  supplied  with  liquid  helium  (He"4). 
Liquid  helium  is  consumed  to  precool  the  He3  returning  to  the  refrigerator. 
Liquid  helium  is  also  consumed  by  the  1.3  Kelvin  coldplate.  Therefore  a 
separate  liquid  helium  container  is  built  onto  one  of  the  copper  shells  of 
the  dewar.  There  is  a thermal  link  (such  as  a thick  copper  braid)  between 
the  cold  plate  and  the  1.3  Kelvin  copper  shield  and  platform. 

The  following  pages  outline  the  thermal  calculations  for  the  1.3  Kelvin 
and  20  mlllikelvln  shields. 
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I.  1.3  K shield  heat  balance 

A.  Heac  flov  into  shield  from  4.2  K CO  1.3  K 

1.  Conduction  chrough  < 10  7 Corr  He  gas; 
q < .7  mW 

2.  Black  body  radiation; 
q .16  mW 

3.  Conduction  through  10  niobium  wires,  4 mils 
(102  microns)  In  diameter  and  300  cm  long; 
q • .05  mW 

4.  Conduction  through  7 support  rods,  fs  inch  (1.27  cm) 
diameter,  1 inch  (2.54  cm)  long; 

nylon  q ■ 1.6  nU 

vespel  SP-22  q ■ 1.4  mW 

B.  Heat'  flov  out  of  shield  due  to  1.3  K cold  plate  of  dilution 
refrigerator;  up  to  20  mW 

C.  Conclusion,  the  excess  cooling  power  of  the  cold  plate  Is  more 
than  enough  to  take  care  of  the  above  sources  of  heat. 

D.  The  thermal  gradient  across  the  copper  shield  from  the  supports 
to  the  cold  plate; 

nylon  AT  - 2.8  mK 

vespel  AT  - 2.4  mK 

II.  20  mK  shield  heat  balance 

A.  Heat  flow  into  shield  from  1.3  K to  20  mK 

1.  Conduction  through  < 10  10  torr  He  gas; 
q < .04  pW 

2.  Black-  body  radiation 
q = .18  VW 

3.  Conduction  through  8 niobium  vires,  4 mils, 

(102  microns)  in  diameter  and  one  meter  long; 

q - . 12  PW 
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4.  Conduction  through  4 vespel  supports, 

.25  Inches  (6.4  mm)  diameter,  2 Inches 
(51  mm)  long;  • 
q - 3.1  yW 

B.  Heat  flow  out  of  20  mK  shifeld  to  dilution  refrigerator  after 
the  system  has  been  cooled;  when  the  refrigerator  Is  at  18  mK, 

Its  cooling  power  Is  3*j  microwatts. 

C.  Net  heat  flow;  the  dilution  refrigerator  can  maintain  the 
experiment  at  20  mK  with  an  expected  heat  leak  of  about  3*j 
microwatts  with  a 2 mK  temperature  drop  along  the  copper 
braid  heat  link.  This  assumes  that  the  heat  generated  by 

, magnetic  levitation  Is  small. 

D.  The  thermal  gradient  across  the  20  mK  shield  is  expected  to 
be  less  than  1 pK. 

III.  Cooldown  time;  the  15  kg  silicon  crystal,  the  copper  shield,  and  the 
ring  suspension  can  be  cooled  from  1.3  K to  20  mK  in  less  than  one 
day. 
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Bar  and  Suspension  Inside  20  mK  S eld 


tZOyr)  Chopper 


LOW  FREQUENCY 
SUSPENSION 


SENSITIVITY  OF  A GRAVITATIONAL  G RADIOMETER  USING  AN  ELASTIC  SOLID 


As  in  most  ocher  measurements,  there  is  a continuous  background  of  signals 
resulting  from  Interaction  of  an  elastic  solid  gradiometer  with  its 

heat  bath,  electronics,  cosmic  radiation,  seismic  noise,  cryogenics 

. Fluctuations  in  this  background  may  resemble  any  kind  of  signal.  Criteria, 
ultimately  based  on  the  theory  of  probability  are  established  to  enable  decisions 

t.  i 

to  be  made,  at  a certain  level  of  confidence,  that  observed  signals  originate 
in  sources  outside  of  the  gradiometer. 

In  Figure  1 is  shown  the  external  signal  which  may  have  any  time 

dependence  whatsoever,  together  will  all  relevant  internal  noise  sources. 

and  are  voltage  and  current  noise  sources  selected  to  account  for  all 

noise  generated  in  the  amplifier  following  terminals  x X . \/,  is  a voltage 

i 

noise  source  describing  most  of  the  thermal  fluctuations  of  the  oy Under  mode, 
is  a current  noise  source  accounting  for  the  noise  of  resistance  /Zt 
which  couples  the  heat  bath  to  output  capacity  and  the  cylinder  mode.  The 
amplifier  is  followed  by  a demodulator  and  low  pass  filter.  The  Maxwell 
network  relations  together  with  integrals  over  the  filter  bandpass  then  lead 

l 

to  expressions  for  the  signal  and  noise  power  outputs,  for  any  given  type  of 
signal. 

Fluctuations 
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The  design  of  the  demodulator  and  low  pass  filter  will  depend  on  the 
dynamical  variables  chosen  to  specify  the  state  of  the  cylinder.  Let  us 


therefore  consider  first  the  available  amplifier  output  voltage.  This  may  be 
written  in  terms  of  the  fluctuating  amplitude  A(tJ/  fluctuating  phase  Cp (“£ ) 
and  mode  angular  frequency  , as 


V(t)  = Aft)  Sit 1 f 


(1) 


We  choose  variables  x and  y defined  by 


X - 


(2) 


*1  * (3) 

to  describe  the  cylinder  state  and  its  time  evolution.  (2)  and  (3)  are  obtained 
from  (1)  by  means  of  electronic  devices  which  multiply  (1)  by  cos  and 

separately  by  Sir\  > fro®  a quartz  oscillator  tuned  as  closely  as  possible 

to  UJ0  . The  separate  channels  then  process  x and  y.  This  multiplication  gives 
for  one  channel 


£Wt/co5  *o£  = /UcnCj?  ■+■  /|  St't,  (2-% t 


(4) 


and  for  the  ocher  channel 


LV(±)sc*  Uj 0t  * A — A CoS  + (£  J 


(5) 


(4)  and  (5)  are  then  coupled  to  low  pass  filters  which  remove  the  terms 
containing  frequencies  £u_>o  and  give  time  averages. 

\ (t’l  ~ ^ l 4 M c-os  (6) 

T ^ = r /1W^  cp[-tU±:  (?) 


r 

The  state  of  the  harmonic  oscillator  representing  the  'cylinder  mode  As 
then  represented  by  a point  in  the  x,  y plane.  This  point  undergoes  a kind 
of  random  walk  as  the  system  evolves. 


Shore  Pulse  Algorithm 


Let  us  assume  that  a-  short  pulse  of  gravitational  gradient  may  produce 
a sudden  displacement  of  the  point  x,  y in  any  direction  whatsoever  with 
equal  probability.  To  compute  the  sensitivity  we  calculate  the  expected  changes 
Ax.  and  4 4^  which  may  be  produced  by  the  background  noise. 

The  power  spectra  of  the  noise  voltage  and  noise  current  _Z"t  are 


given  by 


2 

’r{cwtr-,; 


ikr 


rIL  = 


On  the  far  right  of  (8)  and  (9)  are  the  limits  for  f Let 

Ji-T 

the  power  spectra  of  the  electronics  current  and  voltage  noise  sources  be 

and  P respectively.  These  sources  may  be  somewhat  correlated.  A 
low  pass  Butterworth  filter  is  assumed  to  be  employed  after  demodulation,  with 
upper  cut  off  angular  frequency  UJ^  , averaging  time  ^ and  power  response 

J ^ • It  is  also  assumed  that  £>  and  defined  by  Q s 

*r*  both^^^.  To  good  approximation  the  quantity 

’ ^ My  then  be  c°®puted  as 


<(**>'* 


i 
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The  first  integral  on  the  right  of  (H)  gives  the  wideband  contribution, 
and  the  second  integral  gives  the  narrow  band  contribution.  The  first  tern 
in  the  bracket  of  the  second  integral  gives  the  narrow  band  contribution  of 
the  noise  sources  1/  and  , and  the  second  term  in  the  bracket  gives  the 
narrow  band  contribution  due  to  the  electronics  current  source  . Carrying 
out  the  integrations  gives 


r I FirVl^Cj-) 

8 1 JT 

, I **Tut  + T Pft  ^ _! 

1 O.L,  _ 

In  (12)  we  have  employed  the  relation: 


\jJ  > = J£_ 

“ *ks  JT 

for  a Butterworth  filter.  It  is  illuminating  to  discuss  the  significance  of 
the  terms  in  (12).  The  first  term  tells  us  that  for  the  wideband  noise  the 
fluctuations  in  the  voltage  amplitude  are  essentially  the  same  magnitude  as 
the  noise  voltage  squared  Itself.  The  two  terms  in  the  narrow  band  noise 
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contribution  may  be  written  in  terms  of  (13)  and  the  € folding  time 


£ 


2.0, 


LLJ 

V 


as 


-nkr  2r. 


A 4.  ~n~  Ace  £-/ 

♦S^Q1 


(14) 


The  electronic  noise  raises  the  mode  temperature  from  T to  the  value 


J /l  r -ZrA>  C / I 
4 C£  kT  ) 


are  a factor 


and  the  fluctuations  in  narrow  band  noise  voltage 
smaller  than  the  squared  narrow  band  total  voltage.  It 


<-~Ut 

is  also  clear  from  (14)  that  reduction  of  the  temperature  T must  be  accompanied 


by  improved  electronics. 

For  these  assumptions  an  optimal  filter  is  one  for  which  the  low  pass 
Butterworth  filter  cut  off  is  chosen  so  that  (12)  is  an  extremum.  Setting 

the  derivative  of  (12)  with  respect  to  lUv  equal  to  zero  gives  Lo^  , the 
optimal  bandwidth  as 


Uj>  - 
u 


kT  ( 


I 


II  Pz*  ^1 


1 Ctl  kT 


^~(Prc 

wV*  / 
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(15) 


Let  us  define  a system  "temperature"  by  the  criterion  that  it  is  the  sudden 
change  in  temperature  of  the  cylinder  mode  alone  which  will  give  the  change 
equal  to  the  mean  square  fluctuations  computed  by  (12).  Let  this  temperature 
be  T' . We  have,  then 


¥li‘iz+Ptt  jX^i]  1 


JT 


'V 
UJ  bj 


fc^kr  ) 


(16) 
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solving  for  T'  and  rearranging  terms  then  gives 

(,  rlhs&ttSj)  (1 

^ iC' Lr  1 J 


(17)  then  tells  us  that  the  "effective  mode  temperature"  is  the 
thermodynamic  temperature  times  the  sum  of  a number  of  factors.  The  first 
factor  is  roughly  the  optimal  filter  averaging  time  divided  by  the  mode  € 
folding  time,  the  second  is  roughly  the  ratio  of  the  total  equivalent  squared 
broad  band  noise  current  divided  by  the  total  squared  (narrow  band)  mode  noise 
current.  This  analysis  is  applicable  to  any  gravitational  gradiometer 
described  by  the  noise  sources  and  equivalent  circuit  of  Figure  1. 
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Coupling  and  Quality  Factor 


Expression  (17)  stresses  the  importance  of  a high  quality  factor 
for  an  elastic  solid  sensor.  (17)  also  makes  clear  the  importance  of  having 
the  largest  possible  coupling  (small  L})  together  with  electronics  with  small 
current  and  voltage  noise. 

A 1400  kilogram  aluminum  cylinder,  which  is  carefully  mounted  , may  have 
Q 500,000  at  room  temperature,  with  no  crystals  or  other  attached  instrumen- 
tation. Very  recently  we  installed  a number  of  PZT-8  (lead  zirconate  titanate) 
crystals  on  an  aluminum  cylinder  having  a diameter  close  to  one  meter  and  length 
155  centimeters,  of  mass  3400  kilograms,  at  room  temperature.  The  loaded  Q 
was  150,000,  C2  ■ 3650  x 10”12  farads,  Cj  - 2 x 10~12  farads,  Lj  ■ 5,000  Henries 
fi  ~ 331  0,  r2  ■ 8 x 106  0.  At  lower  temperatures  Professor  J.-P.  Richard  and 
Dr.  Wn.  Davis  obtained  loaded  Q 106  for  an  aluminum  cylinder  with  a parametric 
capacitor.  Dr.  Suzuki,  Dr.  Tsubono  and  Dr,  Hirakawa  have  observed  much  higher 
Q's  with  type  5056  aluminum  at  low  temperatures. 
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